Lifetime statistics of quantum chaos studied by a multiscale analysis 
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In a series of pump and probe experiments, we study the lifetime statistics of a quantum chaotic 
resonator when the number of open channels is greater than one. Our design embeds a stadium 
billiard into a two dimensional photonic crystal realized on a Silicon-on-insulator substrate. We 
calculate resonances through a multiscale procedure that combines graph theory, energy landscape 
analysis and wavelet transforms. Experimental data is found to follow the universal predictions 
arising from random matrix theory with an excellent level of agreement. 
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Introduction. — In a seminal paper of 1917, Albert 
Einstein provided a remarkable insight about the failure 
of energy quantization schemes when applied to noninte- 
grable systems and raised the question "does chaos lurk 
in the smooth, wavelike, quantum world?" [l], This 
problem was practically ignored until 1970, when theo- 
retical physicists led by M. Gut z wilier investigated the 
implications of classical chaos for semiclassical quantum 
systems and founded the research field known as quan- 
tum chaos [3]. Besides quantum mechanics, the problem 
posed by Einstein is of great importance due to the large 
number of systems that can be treated semiclassically 
[H, HI . Among them, electromagnetic waves have stirred 
particular interest. In electrodynamics, quantum chaos 
originates from the isomorphisms between Schrodinger 
and Maxwell's equations in two spatial dimensions and 
is manifested in dielectric optical microresonators whose 
forms mimic classically chaotic billiards. In these geome- 
tries, semiclassical methods and Random Matrix Theory 
(RMT) provide theoretical predictions that have found 
to well agree with several experiments, leading not only 
to fundamental discoveries but also to novel concepts in 
the field of laser devices One of the most funda- 

mental aspects of quantum chaos lies in the universality 
of its eigenmodes, first conjectured by Bohigas, Gian- 
noni and Schmit in 1984 for closed systems [13| and then 
theoretically extended to open media through random 
matrix theory 3, Hj. According to RMT, the spectral 
resonances of a quantum chaotic system show universal 
probability distributions, which depend only on the sym- 
metries of the original equation and the number of open 
channels M considered [ill, QjJ. However, while there 
exists a large experimental literature on the study of the 
position of each resonance (i.e., the real part of the eigen- 
value of an eigenmode), experimental work on the life- 
time statistics (i.e., the imaginary part, or equivalently, 
the resonance width) is still at the beginning. The mea- 
surement of the lifetimes, in fact, is an extremely chal- 
lenging task due to resonances overlapping in the eigen- 
modes spectrum. In this field of research, the present 




FIG. 1. Color Online. SEM image of the open stadium bil- 
liard realized within a two dimensional photonic crystal. 



literature consists only of two experiments, which have 
been performed on microwave billiards characterized by a 
single open channel [HI, [HI . Lifetime statistics extracted 
from measurements have reported to significantly deviate 
from the prediction of RMT for J\f = 1, showing a behav- 
ior whose physical origin is still not understood [18]. Be- 
sides that, work in correspondence of M > 1 (where the 
form of the statistics is strongly different from the case 
AT = 1) is missing, as well as an experimental analysis at 
optical wavelengths, which are extremely interesting due 
to the possibility of exploiting nonlinear effects to control 
quantum chaos phenomena (see, e.g., foL fl2|). 
In this Letter, we study the lifetime statistic of a chaotic 
optical microresonator in the case of M > 1. In a series of 
pump and probe experiments, we injected a broadband 
source signal into a stadium microresonator and collected 
the output by a fiber spectrograph. Resonance lifetimes 
have then been extracted from the power spectrum by a 
multiscale procedure which combines wavelet transforms 
and energy landscape analysis [2o|, HH . We performed an 
experimental campaign on five different samples, and an- 
alyzed more than 700 resonances. Experimentally mea- 
sured statistics were found to perfectly match the uni- 
versal predictions of RMT. 
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FIG. 2. Color Online, (a) Experimental power density spec- 
trum V(X) retrieved for a sample with / = 2.88; (b) enlarged 
portion (a) for A G [1.546/xm, 1.564/im]. 



Sample & experimental setup. — Experiments on 
wave chaos have been traditionally performed in the mi- 
crowave regime due to possibility of using large geome- 
tries (~ mm) and low loss metals acting as mirrors. At 
optical wavelengths, conversely, metals are lossy and low 
absorption dielectrics should be employed. However, a 
naive implementation of a chaotic billiard with dielectrics 
results into large fresnel losses, and RMT statistics would 



no longer be universal [22|. In order to overcome this 



problem ad design an optical cavity that mimics a mi- 
crowave billiard , we resorted to Photonic Crystals (PhC) 
technology 23|, |24[. In particular, we embedded a fully 
chaotic stadium geometry [25] in a two dimensional pho- 
tonic crystal made by a triangular lattice of cylindrical 
air holes in silicon (Fig. [1]). The sample was fabricated 
on a Silicon on Insulator substrate, with a SiC>2 buffer 
layer of 2\im and a Silicon capping layer of 220nm. The 
PhC pattern was written with a standard lithographic 
procedure, involving exposure with an e-beam system, 
development of the photoresist and transferring of the 
pattern on to the substrate with reactive ion etching step 
(see e.g. [26]). The confinement in the transverse direc- 
tion is guaranteed by total internal reflection, whereas 
the PhC cavity ensures light confinement in the plane. 
The sample was designed to work with the electric field 
parallel to the surface (TE modes). Inside a photonic 
crystal, wavelengths A that fall in the PhC gap AA can- 
not propagate due to a zero density of electromagnetic 
states [23|. In such a frequency range, light energy en- 
tering into the billiard gets totally reflected by the PhC 
walls and escapes in the input and output waveguides. 
Other sources of outcoupling losses, besides the waveg- 
uides, may be identified in the Silicon-air and Silicon- 
Si02 interfaces, owing to the finite dielectric permittiv- 
ity of the Silicon, and in the PhC itself. The latter, due 
to its finite extension, is expected to sustain a finite — 
albeit evanescent ly small — energy leakage in the plane. 
All these five elements (input/output waveguides, PhC 
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FIG. 3. Color Online, (a), (b): continuous wavelet transform 
of the power density spectrum V(X) of Fig. [2^ and Fig. EJd, 
respectively. 



and upper/lower interfaces) are completely independent 
and work as J\f = 5 autonomous open channels. 
In order to collect a sufficiently large statistics, we de- 
signed different samples with the same billiard area (of 
700/im 2 ) and diverse filling factor / = ^, being r the 
radius of the single air hole and a the lattice spacing 
(Fig. [1]). By varying the filling factor in the range 
/ G [0.25,0.3], we guaranteed a photonic bandgap for 
electromagnetic TE modes of AA > 400 nm, which is 
sufficiently larger than the bandwidth of the laser source 
employed for the experiments. The experimental setup 
consisted of a polarized C+L band amplified spontaneous 
emission source (central wavelength 1575 nm, bandwidth 
110 nm), which injected light into a single mode fiber that 
was coupled to the sample by a 60x aspheric lens, with 
antireflection coating. Light emerging from the output 
channel was then collimated by a 40x aspheric lens into 
a second single mode fiber, split between a photodetector 
(to monitor the coupling optimization) and an ANDO op- 
tical spectrum analyzer. All spectra were acquired with 
a resolution of 10 pm. 

Multiscale analysis of the resonances. — Figure ((2]) 
shows a typical power density spectrum V(X) obtained 
from a single measurement. We observe strong cluster- 
ing effects, with several resonances overlapping in the 
spectrum (see Fig. [2)3). We emphasize that, due to such 
overlapping, no standard method of resonances extrac- 
tion produced satisfactory results. For this reason, we 
developed an original approach that combines ideas from 
multiresolution analysis and graph theory. More specif- 
ically, we begin by applying a continuous wavelet trans- 
form (W V)(C: s) to the power density spectrum V: 



(wpxc,x)= r v{\)- 

J — 00 V / 



Vx \ x 



(i) 



being \i C scaling and translation parameters, respec- 
tively, and ip the so-called mother wavelet, which is a 
compact function in L 2 (R) possessing a zero mean value 

and satisfying ^ffi df < 00 (^ denotes the Fourier 
transform of ijj). To extract Lorentzian-like linewidths 
out of the spectrum, we choose the following symmet- 
ric wavelet if) = (2 — 4t 2 )e _£ , which is defined from 
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FIG. 4. Color Online, (a) uphill pathways analysis (runners 
as solid lines and nodes as circle markers) of the wavelet trans- 
form of Fig. [3]3, reported for reference as pseudocolor plot; 
(b) theoretical reconstructed spectrum versus experimental 
measure. 



the second derivative of a Gaussian function. The in- 
tegral transform (pQ), when applied to the power density 
spectrum resulting from a measurement, provides a ge- 
ometric visualization of the inner structure of the spec- 
tral resonances — given by the regions of high density in 
Fig. [3] — which change according to the scale x consid- 
ered. As x increases, resonances found at smaller scales 
evolve into tree-like structures and eventually join to- 
gether (Fig. [3}d). The topology of the surface (W V) 
encodes the structure of all the resonances of the sys- 
tem. Isolated resonances, in particular, are observed due 
to the missing of any tree structure in (W V), while a 
link between high-density regions is the signature of the 
presence of a cluster. We describe clusters and isolated 
resonances by an uphill landscape analysis [2l| applied to 
the wavelet transform. Figure 2^ illustrates this analysis 
when applied to the portion of the spectrum of Fig. [2)3. 
In particular, we employ a series of runners (Fig. [4^, solid 
lines), which evolve along x following the curve of min- 
imum steepness on the surface (W V), and nodes (Fig. 
Hk, circle dots) denoting the intersection among different 
runners. We geometrically describe tree-like structures 
among paths of different runners by building an adja- 
cency matrix ^4, with Aij = 1 if the path of runners i and 
j overlaps, and Aij = otherwise. The undirected sim- 
ple graph (5, arising from ^4, carries the information on 
both the structure and the composition of all the clusters 
of resonances in the spectrum. In conjunction with the 
wavelet transform (W V), it provides the required knowl- 
edge to extract position and spectral linewidth of the bil- 
liard eigenmodes. Isolated resonances, characterized by 
= for j ^ z, are described by following the path of 
the corresponding runner up to the first maximum in the 
surface (W V), which yields the highest overlap accord- 
ing to the L 2 distance denned by (pQ). At the maximum 
point in the plane (x, C)? the Full Width Half Maximum 
(FWHM) linewidth SX is ^2 log 2 times the wais t of the 
gaussian part of the wavelet, and reads SX = xy^og2/2 : 
while the resonance position Ao is given by Ao = C I n 
the presence of resonance clusters, which show nonzero 
connections in the adjacency matrix (i.e., A^ ^ for 
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FIG. 5. Color Online, (a)-(b) Position A and FWHM 
linewidth ^A of chaotic resonances extracted by the multi- 
scale analysis, (b) shows a magnified version of (a) in the 
region where the resonances exhibit a random behavior. 



j ^ z), we found optimal (Ao,5A) by a nonlinear sim- 
plex optimization search in the subspace spanned by the 
runners path in the plane (xX) °f the surface (W V). 
Stemming from the theoretical parameters (Ao,5A) ob- 
tained through the combined analysis on (W V) and (5, 
we finally apply a global nonlinear least-square optimiza- 
tion across the entire spectrum for further increasing the 
fidelity of our estimated results. The outcome of this 
procedure are quite satisfactory. Figure HJd compares a 
portion of the theoretical spectrum reconstructed by the 
method described above and the experimental measure. 
As seen, the experimental data is barely visible due to 
an excellent agreement with the fitted curve obtained 
through the multiscale analysis. We highlight that all ex- 
perimental spectra (not shown here), were reconstructed 
with the same degree of precision. 

Discussion of results. — Figure [5] reports the distri- 
bution of the resonances in the plane (Ao, SX) for the five 
samples considered in our experiments. The figure shows 
a clear separation of the resonance dynamics into two 
different regions. A few, isolated short living resonances, 
and many long living modes with positions that are com- 
pletely random in (5 A, Ao). The former, associated to bil- 
liard's bouncing-ball eigenmodes corresponding to nodal 
lines along horizontal and vertical spatial directions, are 
normally ignored in the statistical analysis of the spec- 
trum [27]. We therefore restrict our attention to long- 
living resonances, and study their statistics by applying a 
cutoff rate of SX < 0.6. In order to compare our results to 
the corresponding quantum- mechanical predictions, we 
begin by expressing the scattering matrix S(£) [H, EH, 
which relates the vector of incoming B and outgoing 
A electromagnetic field amplitudes by B = S A, with 
S(£) = 1 - 2iriW(£yV(£y 1 W(£), £ playing the role 
of a quantum-mechanical energy, W the operator model- 
ing the coupling with the environment (i.e., the channel 
space) and V = £ - Ho + zttWW 1 ", being Ho = V 2 + k 2 
the Hermitian Hamiltonian corresponding to the close 
billiard. Resonances of the system originate from the 
poles of the scattering matrix 5, and therefore, from the 
eigenvalues £ = v — i*y of — inWW^ . The eigenvalues 
imaginary part 7, which corresponds to the distribution 
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FIG. 6. Color Online, (histogram) Distribution of resonance 
imaginary part p(j) and theoretical prediction arising from 
RMT (solid line). 



of the resonances width of the open chaotic system, can 
be related to the electromagnetic measured data (Ao, SX) 
from the relation k 2 = ^ = £, and reads 7 = SX/Xq 
(omitting inessential proportionality constants). Proba- 
bility distributions of eigenvalues £ arising from Hamil- 
tonians Ho — inVVW^ have been studied in great detail 
and reviewed by Fyodorov and Sommers in [30|. Closed 
form expression for the probability of decaying constants 
^(7) can be obtained by the supersymmetry method and 
in the limit of weak losses is expressed by the following 
universal distribution: 



P(l) 



ff 



r(.A0 



(2) 



with g = (A/"(7)) _1 . The distribution © is uniquely 
determined by the number of channels A/*, which defines 
the dimensionality of the channel space. When J\f = 1, 
Eq. (|2]) collapses to a simple exponential distribution, 
while for M > 1 it acquires a radically different shape 
characterized by a nonzero probability maximum. Figure 
[6] compares the experimentally measured distribution of 
7 with the theoretical prediction of RMT with M = 5. 
As seen, an excellent level of agreement is found between 
theory and experiments. It is worthwhile remarking that 
we did not employ any fitting parameter in Fig. [6j but 
just the universal distribution (j2j) with M = 5. 

Conclusions. — Motivated by a question posed by 
Einstein in the last century, we investigated the universal 
behavior of quantum systems in the presence of determin- 
istic chaos. By exploiting an isomorphism with Maxwell 
and Schrodinger equations, we designed and experimen- 
tally characterized photonic crystal resonators that are 
fully equivalent to open quantum billiards. Our experi- 
mental results, analyzed within a theoretical framework 
that combines ideas from mult iresolut ion analysis, unam- 
biguously demonstrate the existence of universal statis- 
tics in the lifetime of quantum chaos. Universality is a 
concept of utmost importance in every branch of physics, 
and its implications at the quantum scale are fundamen- 



tal not only to add a piece to the puzzle initiated by 
Einstein's intuition, but also to give rise to new quantum 
devices that exploit universal phenomena and operate re- 
gardless of their microscopic details. 
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